Striped states in weakly trapped ultracold Bose gases with Rashba spin-orbit coupling 
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The striped state of ultracold bosons with Rashba spin-orbit coupling in a homogeneous infinite system 
has, as we show, a constant particle flow, which in a finite-size system would accumulate particles at the 
boundaries; it is thus not a physical steady state of the system. We propose, as a variational ansatz, a 
condensate wave function for a weakly trapped system which behaves similarly to the striped state near 
the center, but does not have particle flow at the boundaries. This state has a line of unquantized coreless 
vortices. We show, by minimizing the total energy, that our modified striped state has lower energy than 
the conventional striped state and it is thus a physically appropriate starting point to analyze striped states 
in finite systems. 



With the recent experimental realization of non-Abelian 
gauge fields in ultracold atomic gases Ol-d], and in par- 
ticular Rashba-Dresselhaus spin-orbit coupling 0, IH, the 
properties of the novel phases possible in these systems 
have come to the fore. Ultracold bosons with Rashba- 
Dresselhaus spin-orbit coupling have been predicted to ex- 
hibit, depending on the details of the interparticle inter- 
action, a plane-wave state, which is a Bose-Einstein con- 
densate (BEC) of a single momentum state, or a striped 
state, which is a BEC made of particles with two opposite 
momenta llq- llOll . At nonzero temperature, exotic phases, 
such as a nematic superfluid phas e lUlll a nd a BEC stabi- 
lized by interactions, are predicted rflaloi . In trapped sys- 
tems, even more exotic phases, such as a half-vortex phase 
and a lattice phase, have been predicted Il4l - [l7ll . When 
the trap potential energy is small compared to the kinetic 
and the interparticle interaction energies, numerical calcu- 
lations predict that in two dimensions phases similar to the 
plane-wave and the striped phases in a homogeneous sys- 
tem may emerge ifTsLITql . 

In a homogeneous system, the striped phase has, as we 
show, alternating local particle current along the stripes, 
which accumulate particles at the boundaries, and thus it 
cannot be a physically realizable steady state in a finite sys- 
tem. In this paper, we address the question of how such a 
striped state can be modified to embed it into a finite sys- 
tem. Two possible remedies are surface currents at the ends 
of the stripes which recirculate the particle flow, or flow 
between the stripes in bulk. We explore this latter solution 
in this paper, considering a system trapped by a weak har- 
monic potential in the direction (along y) of the stripes. For 
simplicity, in this first step, we do not confine the system 
transverse to the stripes. We propose a variational ansatz 
for a modified striped condensate wavefunction which be- 
haves as the striped state near the center, but does not have 
a net current flow at the boundary. The wavefunction is 
Gaussian in the direction of the stripes, with a complex 
phase describing divergence-free particle flow. Fixing the 
spatial extent of the Gaussian by minimizing the energy, we 
show that the modified striped state has lower energy than 
a simple striped state, and thus is a physically acceptable 



candidate for the ground state. 

Local current in the striped phase. We consider bosons 
in two hyperfine states, labeled f and \., with pure Rashba 
spin-orbit coupling, described by the Hamiltonian 
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where ^(r) = (^(r), if)±(r)) is the particle annihilation 
operator, V(r) is the potential of the confining trap, and 
A = — k(o~ x , a y , 0) is the non-Abelian vector potential, 
with coupling strength k. We assume s-wave interactions 
between particles, 
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The Hamiltonian without interactions has a circle of 
degenerate single-particle ground states with momenta 

\jPx + Py = K anc ^ Pz = 0- The striped state, a super- 
position of two states with opposite momenta on the circle 
of degenerate single-particle ground states, is the preferred 
ground state in the homogeneous (V(r) = 0) infinite sys- 
tem, if we take the renormalization of the interaction into 
account Il0lll8lll9ll . The macroscopic wavefunction of the 
striped state with momenta (±k, 0, 0) is 
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-i sin(Kx) 
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where n is the density of particles. 

The particle current in the presence of the gauge field, 
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does not vanish; there is a local flow of particles in the 
striped state all the way out to y = ±oo. (We set h = 1 
throughout.) Figure [TJ plots (j x (r),j y (r)) as a function 



2 



of kx and ny. The flow in the y direction alternates as 
we move along the x axis. In a finite-size box, this flow 
would drive particles to the boundaries in the y direction, 
accumulating particles at the edge. Thus the striped state 
per se is not a stable steady state in a finite-size system. On 
the other hand, experimentally realizable ultracold atomic 
systems are inevitably finite, and we are in need of finding 
a physically acceptable state which is similar to the striped 
state. 
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FIG. 1: The particle current j(r) of a striped state as a function 
of x and y in units of k . The thickness of the arrows indicates 
the magnitude of the current. 



The particle circulation in the striped state is not quan- 
tized, since the local velocity, v st = j st /n is 

/ 2k 2 \ 

Vxv„= 0,0, cos(2kx) . (5) 

\ m J 

Such nonquantized vorticity is a generic feature of systems 
coupled to external gauge fields, e.g., type-II superconduc- 
tors. 

Nonvanishing local flow is a common feature of striped 
states in general Rashba-Dresselhaus spin-orbit coupling 
fields of the form A = — k(o~ x , r\o v , 0), where < rj < 1. 
When rj = 0, i.e., an equal mixture of the Rashba and Dres- 
selhaus spin-orbit couplings, the local current vanishes. In 
the experimentally realized spin-orbit coupling to date is 
an equal mixture of the Rashba and Dresselhaus spin-orbit 
couplings [1-3] (with no problem of particles accumulat- 
ing at the boundaries). 

Modified striped state. We now assume that the sys- 
tem is confined in the y direction by a harmonic potential, 
V(r) = muj 2 y 2 /2. Taking the striped state © as a start- 
ing point, we assume a macroscopic wavefunction for a 
modified striped state of the form 



* ms (r) = 7 (x,y)e < *^ ( cos ^ n 



(6) 



where j(x, y) and (j)(x, y) are real functions of x and y. 
Going beyond this form of ansatz, for example by modify- 
ing the spin part of the wavefunction, is beyond the scope 
of this paper. In this state, the particle density is j 2 , and 
the particle current is 



j(r) = — [V</> + (0,-Ksin(2/cx),0)] 
m 
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Two physical conditions on ^ ms are that the current van- 
ishes at large \y\ and, in equilibrium, V • j = 0. A natural 
choice for a system harmonically confined in the y direc- 



tion is 
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where £ is a positive dimensionless variational parameter, 
and the normalization c is determined by fixing the total 
number of particles N, 
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Here L is the linear size of the system in the x and z di- 
rections, and n ms = j dyj 4 / J dyj 2 is a mean particle 
density. With this form for 7, the condition V • j = is 
satisfied for 
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Then, 



j(r) 



K 2 7 2 



{Qny cos(2kx), — sin(2Kx), 0) . (11) 
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In the limit Q — > 0, we recover the striped state (0 and 
current ©. Figure |2]plots the current for £ = 1 as a func- 
tion of kx and ny. When ny is small, the current behaves 
similarly to that in the striped state and flows in the ±y di- 
rections. However, as \y\ grows, the current bends towards 
the ±x directions, and particles no longer accumulate at 
the boundaries in the y direction. Near the center, this state 
in a finite system behaves similarly to a striped state. 

The modified striped state has a line of vortices, along 
the z-axis, at r = (utt/2k, 0, 0), where v takes integer 
values. While the particle velocity at the center of the vor- 
tices vanishes, the density remains nonzero; the vortices 
are coreless. The curl of the local velocity in the modi- 
fied striped phase, v = j/7 2 , is simply that of the original 
striped phase, given by Eq. d5j [2011 . 

Energy. We turn now to examine when the modified 
striped state is energetically favorable at zero temperature. 
In the mean field, the energy per particle of the modified 
striped state is 
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FIG. 2: The particle current j(r) of the modified striped state 
©, with ( = 1, as a function of a; and y in units of k . The 
magnitude of the current, shown by the thickness of the arrows, 
falls off at large \y\. Note the line of vortices along the x-axis. 



where = k 2 /2m, 
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and e t = raw 2 /4k 2 . 

First we compare the kinetic and interaction terms in 
braces in Eq. (fT2ll with the corresponding terms of the sim- 
ple striped state (0, whose energy per particle is 
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(E) s = — (3g n + 3g u + 2g n ) 
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The kinetic energy of the modified striped state is positive 
and thus larger than that (zero) of the simple striped state. 
Furthermore, the ratio of the interaction energy in the mod- 
ified striped phase to that in the simple striped phase is 
n ms j n. More relevant is to compare the energy of the mod- 
ified state with a straightforward extension of the striped 
state with a Gaussian-striped state described by the con- 
densate wavefunction, 
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Although this state is a natural extension of the striped state 
with a Gaussian decay, the divergence of the current does 
not vanish, and is not physically acceptable. Nonetheless, 
since the correct ground state must have lower energy than 
the Gaussian-striped state, we compare the energies of the 
Gaussian-striped and the modified striped states. The en- 
ergy of the Gaussian-striped state is 
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(16) 



The energy difference between the modified striped state 
and the Gaussian-striped state arises just from the kinetic 
energies, 



(E) ms — (E) 
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Thus, modifying the striped state with a spatially depen- 
dent phase to make the current divergence-free decreases 
the energy below that of the Gaussian extension of the 
striped state, showing that our ansatz wavefunction for 
the modified striped state is a promising weakly confined 
striped-like ground state. 

We also compare the energy of the modified striped state 
with that of a Gaussian extension of the plane-wave state, 
described by the wavefunction 



* = ce -c»v/2f_ ( 1 , 
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(18) 



The current vanishes in this phase and it is thus physical. 
Its energy is 

\ E )gp = £ fe 2 + 8 ^ L2 Wt + 9u + 2 9n) + ^ • 

(19) 

The kinetic energy and trap energy are the same as for the 
Gaussian-striped state. The interaction energy of the plane- 
wave state is lower when g^ < (g^ + g^)/2, in which 
case the mean-field ground state in an infinite system is also 
a plane- wave state [8]. Thus, within mean-field coupling 
with weak confinement, the modified stripe phase is pre- 
ferred over the Gaussian-plane-wave ground state if (and 
only if) the interspecies interaction is larger than the aver- 
age intraspecies interaction. 

Finally, we determine £ variationally; the minimization 
of {E) ms with respect to £ implies that the optimal £ is a 
function of the dimensionless ratios e k /e t and €i/e t . Fig- 
ure [3] shows the resulting energy of the modified striped 




FIG. 3: The energy of the modified striped state (E) m s as a func- 
tion of the kinetic and interaction energies; all energies are in 
units of et . 
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FIG. 4: The inverse width, f, of the Gaussian 7, as a function 
of the kinetic and interaction energies in units of et; note that 
the figure is turned 90° from Figs. [3] and [5] The surface rises 
indefinitely as the kinetic and interaction energies become small. 



state, {E) ms /e t , as a function of t^je t and 6j/e t , while 
Fig.|4]shows the optimal £. 

As to — > 0, the energy (fT2l l is minimized when £ = 0, 
which is the simple striped state; as can be seen in Fig. [4] 
£ grows with increasing u;. In the limits of kinetic energy 
dominant over the interaction energy, and vice versa, we 
can obtain the asymptotic behaviors of £ and the corre- 
sponding energy analytically. 

First, when e k 3> e t , by minimizing (fT2] | and ignor- 
ing the interaction term, one obtains ( ~ 2y/ e t /ck = 
\^2moj/K 2 . The energy to leading order is (E) ms ~ 
uj /2y/2. This asymptotic behavior is valid when the es- 
timated interaction energy using the asymptotic value of 
£ is smaller than the kinetic and the trap energies, which 
yields the condition ei/e t <C (efc/e t ) 3 / 4 /2. In the other 
limit, 6i 3> a similar analysis leads to £ ~ (2e t /ei) 2 ^ 3 
and (E) ms ~ 3(e^e t /4) 1/3 . The condition at which this 
asymptotic behavior is valid is ej/e t 3> (ek/e t ) 3 ^ 4 /3. In 
this limit, the extent of the Gaussian in y is independent 
of k. For typical experimental parameters of spin-orbit 
coupled 87 Rb [2], u/2ir ~ 50 Hz, k ~ \/2vr/800nm, 
/V ~ 5 x 10 5 , and L ~ 100/wn, one has e fc /e t ~ 10000 
and €i/e t ~ 2300; here one is in the latter limit, with 
C ~ 0.008. 

The optimal value of £ for the Gaussian-striped phase 
is similarly a function of the dimensionless ratios e k je t 
and £j/e t . In Fig. \5\ we show the difference of the en- 
ergy of the modified striped state and the Gaussian-striped 
state in a trap, as a function of ek/e t and ei/et- The energy 
of the modified striped state always lies below that of the 
Gaussian-striped state, with the difference becoming larger 
as the kinetic energy becomes large compared to the trap 




FIG. 5: The energy difference of the modified striped state and 
the Gaussian-striped state, (E) ms — (E) gs ; again all energies are 
in units of et . 



energy. 

In summary, the modified striped state is not only physi- 
cally allowed, but it is also an energetically favorable start- 
ing description of steady-state striped states in a harmonic 
trap. The arguments given in this paper hold in two as well 
as three dimensions, mutatis mutandis, and thus the state 
© is the preferred striped-like ground state in either situ- 
ation. As we noted earlier, the modified striped state has 
a line of coreless vortices, which can be probed to identify 
the state in possible future experiments. 

We are grateful to Jason Ho for illuminating discussions, 
and to Xuebing Luo for pointing out a simple error in 
Eq. (12) of the published version of this paper. This re- 
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